The effects of superconductor-stabilizer interfacial resistance on quench of a pancake 

coil made out of coated conductor 

G. A. Levin, W. A. Jones and P. N. Barnes 

Air Force Research Laboratory, Propulsion Directorate, 

Wright-Patterson Air Force Base, OH 45433 USA 



O 

(N 



oo 

(N 

O 

o 
u 

CO 



T3 

C 

o 

(N 
> 

OO 

en 

i> 

o 
o 



X 



K. A. Novak 

Department of Mathematics, Air Force Institute of Technology, 

Wright-Patterson Air Force Base, OH 45433 USA 

(Dated: July 30, 2010) 

We present the results of numerical analysis of normal zone propagation in a stack of 
Y Ba^CuzOi - x coated conductors which imitates a pancake coil. Our main purpose is to deter- 
mine whether the quench protection quality of such coils can be substantially improved by increased 
contact resistance between the superconducting film and the stabilizer. We show that with increased 
contact resistance the speed of normal zone propagation increases, the detection of a normal zone 
inside the coil becomes possible earlier, when the peak temperature inside the normal zone is lower, 
and stability margins shrink. Thus, increasing contact resistance may become a viable option for 
improving the prospects of coated conductors for high T c magnets applications. 

PACS numbers: 84.71. Ba, 84.71.Mn, 74.72.-h 



I. INTRODUCTION 

High temperature superconducting wires such as 
YBa^Cu'iO'j-x (YBC'O) coated conductors [l[ have a 
potential for a number of important applications. All of 
them, and especially the high field magnets, require ef- 
fective quench detection and protection systems 2]. Cur- 
rently manufactured coated conductors are characterized 
by very low speed of normal zone propagation (NZP) 3- 
5]. This makes it difficult to promptly detect the onset 
of a quench with the result that in a poorly cooled (adia- 
batic) coil the temperature of the hot spot will rise above 
the safe limit before the quench protection system is en- 
gaged, irreversibly damaging the conductor. Unlike the 
conventional low— T c wires in which a quench can be eas- 
ily triggered by a small amount of energy, the coated con- 
ductors are substantially more stable. Thus, the quench 
in devices made out of coated conductors has a potential 
to become a Black Swan event - a rare, unpredictable 
occurrence with catastrophic consequences. 

A promising approach to improving the quench pro- 
tection quality of coated conductors is to increase the 
contact resistance between the superconducting film and 
stabilizer. In Refs. [y, 0] it was shown, using a model 
of a linear conductor, that the NZP speed can be sig- 
nificantly increased by placing a high resistance interface 
between the YBCO film and copper stabilizer. The effect 
of increased contact resistance on the superconducting 
wire has been studied in the past. A practical reason for 
experimental and theoretical work in this area was a po- 
tential for reduction of coupling AC losses in the low T c 
wires. A strong decline in stability margins which results 
from large contact resistance and technological difficul- 
ties of making such resistance consistently uniform pre- 
vented the incorporation of this approach into practice. 
Earlier patents on the subject are difficult to locate. An 



extensive, although perhaps not entirely comprehensive, 
list of references can be found in the reviews [g, Q • 

The property of high contact resistance to increase the 
NZP speed was not addressed in this earlier research 
probably because increasing the NZP speed in low T c 
wires was not necessary as it was sufficiently high to as- 
sure adequate quench protection. The advent of coated 
conductors has changed the priorities. Coated conduc- 
tors are much more stable than the low T c wires, so that 
the reduction of stability margins accompanied by in- 
creasing NZP speed and with it an earlier detection of 
quench may be desirable. The sequential stages of the 
manufacturing process of coated conductors are also con- 
ducive to incorporating an additional step of adding a 
thin, high resistance interface between the YBCO film 
and stabilizer. 

Here we present an analysis of a model of quench in a 
pancake coil made out of coated conductor. It should be 
noted that increasing NZP speed is not in and of itself 
our goal. The propagating quench builds up the electric 
potential difference across the coil. Typically a quench 
protection system is engaged when this potential exceeds 
a certain threshold. Therefore, a goal of a quench detec- 
tion system is to reliably detect quench in the shortest 
amount of time, thus minimizing the peak temperature 
inside the normal zone. The effect of increased contact 
resistance on quench is not straightforward. The NZP 
speed increases, but the power dissipation at the front of 
the propagating NZ also increases. Thus, in the worst 
case scenario of adiabatic quench these two effects com- 
pete with each other. We show that the increased contact 
resistance has an overall positive effect by increasing the 
rate at which the potential drop across the coil increases 
and, at the same time, reducing the rate at which the 
peak temperature inside the normal zone rises with the 
potential drop across the coil. Thus, at the moment of 
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FIG. 1: (a) A sketch of a coated conductor (not to scale). 
A thin superconducting (YBCO ) film (3) is deposited on 
metal substrate (2). The resistive interface (4) segregates cop- 
per stabilizer (1) and YBCO film, (b) A stack of TV coated 
conductors (each similar to that in (a)) segregated by elec- 
tric insulation with finite thermal conductivity. The stack as 
a whole is thermally insulated from the environment (adia- 
batic). A blob in the center of the stack indicates a nucleated 
normal zone, (c) A model described by Eq. (2) corresponds 
to all conductors in the stack being connected in series. The 
arrows indicate the lateral heat transfer (along the conduc- 
tor), as well as the thermal coupling over the finite distance 
2L which corresponds to the heat exchange in the transverse 
direction, between the neighboring conductors in the stack 
shown in (b). 



quench detection the peak temperature inside the coil 
will be lower, which leaves a greater safety margins for 
the quench protection system. 



II. MODEL OF QUENCH PROPAGATION IN A 
PANCAKE COIL 

In [1(| it was shown that the 3D equation of heat con- 
duction in a thin tape-like composite wire, Fig. 1(a), can 
be reduced to a 2D (planar) equation if the heat flux from 
the surface does not greatly exceed 1 W/cm 2 . Then, the 
temperature variation across the thickness of the tape 
constitutes a fraction of a degree, much smaller than the 
variation of temperature along the wire. The ID approx- 
imation is valid as long as the thermal diffusion length is 
greater or comparable to the conductor width. In coated 
conductors this condition is met. 

A quench in pancake coil can be adequately modeled 



by a stack of coated conductors segregated by layers of 
insulation as shown in Fig. 1(b). One can view such a 
stack as part of the straight section of a race-track pan- 
cake coil. Incidentally, the stacks of coated conductors 
were used in experimental studies of quench reported in 
[5| . The model of the coated conductor itself is the same 
as in Ref. Q. We take into account four constituencies 
of the coated conductor: stabilizer, substrate, supercon- 
ducting YBCO film, and infinitesimally thin interface be- 
tween YBCO and stabilizer which accounts for the finite 
contact resistance, Fig. 1(a). 

The system of equations describing the heat exchange 
among the coated conductors in the stack takes form: 



C 



m 
at 



d_ 

dx 



K 



.d_T\ 

dx 



= Qi-K ± (Ti(x)-T i+1 (x)) (1) 

-IdCZHaO-IUifc)). 

The conductors in the stack, each of length 2L, are enu- 
merated from i = 1 to N. The dimensionality of Eq. (JTJ 
is [W/cm 2 ]. Here K± = k ins /2d ins [Wcm^ 2 K^ 1 ] is 
the heat conductance across the insulation layer seg- 
regating the neighboring turns of the coil. We con- 
sider that both sides of the conductors are wrapped in 
an insulating film of thickness di ns . Equation (1) de- 
scribes the worst case scenario, from the quench pro- 
tection point of view, of fully adiabatic case with no 
heat exchange between the coil and the environment. 
C = c\di + c 2 d 2 + 2ci ns d ins [Jcm~ 2 K^ 1 ] is the com- 
bined heat capacity of the coated conductor plus insu- 
lation per unit area, K = k\d\ + k 2 d 2 [W K~ x ] is the 
effective in-plane thermal conductance (the in-plane ther- 
mal conductance of the insulation can be neglected). The 
subscripts 1, 2, and ins refer to the stabilizer, substrate, 
and insulation, respectively. The thicknesses of the stabi- 
lizer and substrate are denoted as d± and d 2 respectively. 
Here Cy,c%, Ci nsi ki,k 2 , and ki ns are the respective mate- 
rial constants. Qi = J_\ qi(z)dz is the density of the 
internal heat sources integrated over the thickness of the 
conductor. 

Since all conductors in the stack shown in Fig. 1(b) 
are parts of one conductor and carry the same current, 
it is convenient to reformulate the problem of N coupled 
equations as one equation similar to that analyzed in 
Ref.Q. We will consider the conductors in the stack to 
be connected in series, as shown in Fig. 1(c). The total 
length of the resulting conductor is 2NL. The difference 
between the quench propagation in a single linear con- 
ductor and that in a coil is the thermal coupling between 
the turns which in our case translates into thermal cou- 
pling between the points separated by a finite distance 
2L, as shown in Fig. 1(c). Thus, the system of Eqs. ([1]) 
can be rewritten as 
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dT 
~dt 



d_ 
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K 



dT 

dx 



Q-K ± {T(x)-T(x + 2L)) (2) 

-K±(T(x)-T(x-2L)). 

Any pancake coil, round or race track shape, can be mod- 
eled by a similar equation, except that the finite distances 



between the coupled points will vary along the wire. 

The second equation defined for a conductor shown in 
Fig. 1(c) describes the instantaneous redistribution of 
current between the superconducting film and stabilizer 
and is determined by the condition of charge conservation 



dx 



Vi-Vs 
R 



(3) 



where Ji is the linear density of current {A cm~ l ) flowing 
through the stabilizer. Vi and V s are the local electric 
potentials of the stabilizer and superconductor, respec- 
tively and R [51 cm 2 ] is the resistance of the unit area 
of the interface (contact resistance). This condition can 
also be used in the form 



ox \ ox 



= E\ — E s . 



(4) 



Here E\ = — dVi/dx and E s = —dV s /dx are the electric 
fields in the stabilizer and superconductor, respectively, 
and we do not assume that R is uniform. 

The integrated area density of heat sources in @ takes 
the form 



Q = ^Ef 

Pi 



R 



J S E S , 



(5) 



where J s = J — J\ is the density of current in the super- 
conductor and J = const is the total transport current 
density. 

The constituent relationship between current density 
and electric field for a superconductor will be used in the 
same form as in [7] 



E„ 



p±( 7 - 

d s y s 


-Jo), 


if J s > J c 


o. 




if J s < J c 



(6) 



where p s is the normal state resistivity and d s the thick- 
ness of the superconducting film. Hereafter, we will adopt 
a linear temperature dependence of J C (T) [11] : 



J c {T) = a{T c ~T); T < T c . 



(7) 



For the stabilizer, the conventional Ohmic relationship 
will suffice at all temperatures: 



El = £ a. 

rfi 



(8) 



In the normal state the resistance of YBCO film is much 
greater than that of the stabilizer, 



d s di 



(9) 



The final step in formulating this model is to 
present (0 and (d]) in a proper dimensionless form. The 
current sharing temperature T% is defined by the condi- 
tion J C {T\) — J. Let us introduce a dimensionless tem- 
perature 9: 



e 



T-Ti 

T c - Ti ■ 



(10) 



Then, Eq. takes the form 

J c = J(l-0). 



(11) 



Let us introduce a fraction of the total current u that 
flows through the stabilizer 



Ji = Ju; J s = J(l - it); < u < 1. 
For 6 < 1 equation © takes the form: 



Es = 



(p s /d s )J(8-u), u<6 
0, u> 6 



For 6» > 1 



E„ 



z J ^- v)i 



Then, Eq. (J3]) can be written in a compact form 

d ( du\ 

— A 2 — = u — r max [0, min(0, 1) - u] , 
ox V ox I 



where 



and 



r = l± -r >i 
Pid s 



X 



Rdx 

Pi 



1/2 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



is the current transfer length which determines the length 
scale of the current exchange between the superconductor 
and stabilizer. 

Taking into account ([3]), the heat sources Q defined by 
Eq. ([5]) take form 



Pi J 2 



<r /?J'M^) +-/"(l-«) 



—u-— (r— 
d\ dx \ dx 

' (18) 

We will measure the distances in units of thermal dif- 
fusion length It and time in units of 7 -1 , where 



h = {DtH) 1/2 \ 7 = PiJ 2 /d!CAT. 



(19) 



Here Dt = K/C is the effective in-plane thermal diffusiv- 
ity of the conductor, AT = T c — Ti, and the increment 7 
determines the characteristic time required for the Joule 
heat generated in the stabilizer to warm the conductor 
by temperature AT. 

In dimensionless variables, Eqs. @ and (|15p take the 
form 



86 d 2 6 _ (du 

dT~~de~ u+r \dii 



(1 U) dd{ r dU (21.) 



-/ 5± (20(0-0(e-2O-0(£ + 2*)); 

d ( du\ 

~8f \ r ~8f I =U ~ T max i ' min (o' ) 1) - u \ ■ 



(21) 



with r = "ft, £ 
Here 



x/l T (-tN < £ < £/V), and •£ = L/h 



Kj_ATdi 

Pi J 2 ' 



(22) 



Notice that Eq. (|2"0")) does not depend on the concrete 
form of the constituent relationship between electric field 
and current density in the superconductor. The specifics 
of the constituent relationship enters only in the charge 
conservation condition given by (UJ) and its dimensionless 
versions ([15)) and (|2ip . 

The relative role of the interface resistance is deter- 
mined by the parameter 



A 2 R_ R = pit. 
Ro di 






KAT 



(23) 



The results will not depend on the value of T as long as 
r> 1. Hereafter, for the purpose of numerical solution, 
we take T = 10 2 . 

The main purpose of our analysis is to determine how 
the peak temperature inside the coil changes with the 
electric potential drop across the coil which is typically 
used to detect the normal zone. As a rule of thumb, the 
length of the normal zone, however defined, is considered 
to be a good proxy for the potential drop. Hereafter we 
will use a more accurate definition of the potential drop 
V across the coil due to developing quench: 

f NL n J P Nl 

v ( t ) = ^Ji{x,t)dx^^-l T u(e,t)d£- 

J-nl «i »i J-Ne 

For simplicity, in this analysis we will neglect the tem- 
perature dependence of the resistivity. The characteristic 
voltage 



V 



PiJ, 



(25) 



determines the potential drop due to NZ with the length 
of the order of the diffusion length. 



III. SPECIFIC EXAMPLE: MATERIAL 
PROPERTIES OF COATED CONDUCTORS AT 

4.2 K 

The solutions of Eqs. (20) and (21) provide a universal 
picture of quench propagation, independent of the values 
of the material parameters. To translate these results to 
a case of concrete materials and operating temperature 
we need to calculate the characteristic constants, such as 
dimensionless parameters r and k± as well as the incre- 
ment 7 and the diffusion length It- One of the promising 
near future applications of coated conductors is magnet 
coils operating at or near liquid helium temperature[3j. 
As an example we will use the following values of the 



material constants at or near 4 K [ll| : 
Copper stabilizer: 

pi « 1.5 x 10~ 8 £7 cm; 

fci « 6.3 W/cmK; 

Cl w 8.1 x 1CT 4 J/cm 3 K; 

d\ = 40/im. 



(26) 



Substrate (Hastelloy): 



P-2 



123 x 10" 6 O cm: 



(27) 



k 2 « 8 x 1CT 4 W/cm K; (k 2 = LT/p 2 ); 

c 2 « 9 x 1CT 3 J/cm 3 K; 

d 2 = 50 pm. 

The substrate thermal conductivity was estimated from 
the value of resistivity using Wiedemann-Franz rela- 
tionship (shown in parenthesis) with L = 2.45 x 
10~ 8 WQ/K 2 . 

Insulation (Polyimide Kapton™): 



C-ins 

"ins 



1.1 x 10~ 3 J/cm 3 K; 

1.1 x 10~ 4 W/cmK; 

di n3 ~ 25^m. 



(28) 



The effective thermal in-plane conductance and the 
heat capacity are given by 

C = cidi + c 2 d 2 + 2c ins d ms w 5.3 x 10~ 5 J J cm 2 K{29) 
K = kidi + k 2 d 2 w 2.5 x 10" 2 W/K. 

It is noteworthy that at liquid helium temperature the 
main contribution to the in-plane heat conduction is pro- 
vided by copper, but the main contribution to the con- 
ductor heat capacity is provided by the substrate. This 
creates a possibility of a substantial temperature differ- 
ence between the stabilizer and substrate at the front of 
the propagating NZ. The model here is based on an as- 
sumption that the heat conduction of the buffer between 
YBCO and the substrate is high enough, so that such 
temperature difference between the substrate and stabi- 
lizer can be neglected. The in-plane thermal diffusivity 



D T = K/C w 470 cm 2 /s. 



(30) 



For the sake of estimate we take the critical cur- 
rent density at the operating temperature To equal to 
J c (Tq) — 300 A cm' 1 and the transport current density 
J = 150 A/cm. Then the current sharing temperature 
Ti = (T c + T )/2. Hereafter we take T = 4.2 K and 
T c « 90 K, so that T x « 47 K and AT = T c -T x = 47 K. 
Then, the increment 7 (Eq. (19)) is equal to 



7 



34 s" 



while the thermal diffusion length 



It = (DtM) 1 ' 2 



KdxAT 
Pi J 2 



3.7 cm. 



(31) 



(32) 



It should be noted that the estimates of It and 7 based 
on the assumption that the material parameters are tem- 
perature independent with the values equal to those at 
T = A K tend to produce a significant quantitative er- 
ror. The resistivity and thermal conductivity of copper 
vary strongly within the superconducting range of tem- 
peratures and beyond it. Our purpose here is simply to 
illustrate the effect of the contact resistance keeping other 
material parameters unchanged. The heat conductance 
across the insulation k± can be estimated as 



A" 



h'; 



2di, 



2.2 x 1CT 2 Wcm~ 2 K~ l . (33) 



Correspondingly, the dimensionless coupling constant 
(Eq. (22)) is given by 

Kx « 12. (34) 

The large value of the dimensionless heat transfer con- 
stant k± is noteworthy. It indicates very effective heat 
redistribution between the turns of the coil. It should be 
noted that in the case of an individual linear conductor 
cooled from the surface the condition of cryostability is 
k± > k c < 1 [lOJ . In such a conductor the condition of 
cryostability means that the heat generated in the stabi- 
lizer cannot maintain its temperature above the critical 
temperature. In an overall adiabatic coil the large value 
of kj_ indicates rapid heat propagation across the insu- 
lation. This large thermal coupling between the turns of 
the coil is due to low resistance of the stabilizer, rather 
than a particularly high thermal conductivity of the in- 
sulation. The Joule heat generated in the stabilizer is 
small so that even relatively small temperature differ- 
ence across the insulation ST w AT/n± is sufficient to 
transfer a large part of the generated heat flux to the 
neighboring sections. Also should be noted that it is 
the electric insulation that provides the main thermal 
resistance to the heat exchange between the turns, not 
the substrate. Indeed, k 2 /d 2 ~ 16 x 1CP 2 WcmT 2 K~ x , 
while K x ~ 2.2 x 1CT 2 Wcra' 2 K~ x (Eq. (33)). This 
is the reason why our approximation of uniform across 
the thickness of the conductor temperature is well justi- 
fied. Below we will consider two different values of the 
heat transfer constant. One is given by Eq. (34) and the 
other, much smaller, kj_ — 0.2. This value corresponds 
to greater thickness of the insulator and/or greater resis- 
tance of the stabilizer. 

The characteristic contact resistance defined by Eq. 
(23) is given by 

#o~ 5.1 x lCT 5 f2 cm 2 . 



(35) 



(36) 



The characteristic voltage, Eq. (25), 
V « 2 mV. 

IV. RESULTS 



The system of equations (|2U)l and (|21l) were solved 
numerically by using an IMEX Crank-Nicolson/Adams- 



Bashforth method in conjunction with a fixed-point 
method to solve the Poisson equation with a nonlinear 
source term. The purpose of the numerical solutions is 
to determine how the voltage drop across the coil ( Eq. 
([24]) ) changes with time for different values of the con- 
tact resistance, how the peak temperature inside the NZ 
changes with time, and, by combining these results, how 
the peak temperature inside the coil changes with the 
voltage drop across the coil. 

The solutions 0(£, t) presented below correspond to the 
initial condition in the form of a Gaussian in the center 
of one of the conductors in the middle of the stack. 



6(Z,0) = (a-6 )e-t 2 / 2s2 +t 



(37) 



In all examples shown below we keep S = y2- The peak 
temperature of the initial perturbation given by the con- 
stant a determines whether the nucleus of the normal 
zone will trigger its expansion, or the extra heat associ- 
ated with the initial perturbation will dissipate through- 
out the stack without triggering NZP. Since we choose, 
as an example, the transport current equal to 50% of the 
critical current, J c (Tq) = 2 J, the dimensionless operating 
temperature, see Eq. (11), 



% = l- J -^ = -l. 



(38) 



For (f2~Tj) we used the following boundary conditions: 



i(£N) = u{-£N) = 0. 



(39) 



The results reported below do not depend on the exact 
form of the boundary conditions because they are ob- 
tained from the solutions of Eqs. (|2T)|) and (|2"TT) for which 
the regions of elevated temperature do not yet extend to 
the boundaries of the stack shown in Fig. 1. 



A. Insulation with low heat conductance («j_ = 0.2). 

Figure 2 illustrates the differences between the growth 
of the normal zone in two stacks of coated conductors. 
Figure 2(a) shows the initial condition, Eq. (37), a tem- 
perature spike in the center of the stack. The stripes, 
both white and gray, indicate the individual conductors. 
The alternating coloring of the stripes allows for a better 
visual resolution. Insulation between the conductors is 
not shown. The vertical axis is the local temperature of 
a conductor. Figure 2(b) shows the temperature distri- 
bution in the stack composed of coated conductors with 
low contact resistance, R/Ro = 0.2, characteristic of the 
currently manufactured coated conductors [7J. Figure 
2(c) shows the temperature distribution after the same 
elapsed time in a stack composed of coated conductors 
with much higher contact resistance R/Ro = 10. All 
other parameters of the model remain the same. No- 
tice the substantial temperature differences between the 
neighboring conductors. This is the result of relatively 
weak thermal coupling kx — 0.2. 




FIG. 2: Temperature profiles (vertical axis) in the stacks of coated conductors obtained from the solutions of Eqs. (20) and 
(21). Each stripe in the figure, of either coloring, corresponds to an individual coated conductor like the ones shown in Fig. 
1(b). The insulation layers between the conductors are not shown, (a) The initial condition - a temperature perturbation in 
the center of the stack that gives rise to NZP. (b) The temperature profile after certain elapsed time in the stack composed 
of conductors with low contact resistance, (c) The temperature profile after the same elapsed time in the stack composed of 
conductors with higher contact resistance. The high contact resistance leads to faster expansion of the NZ. It also leads to 
higher peak temperature. In both cases the thermal coupling is characterized by k± = 0.2. 



These results illustrate the apparent inherent contra- 
diction laying at the foundation of the concept of high 
contact resistance. In the stack made out of coated con- 
ductors with high contact resistance the normal zone 
propagates manifestly faster, both in lateral and trans- 
verse directions. However, the amount of heat generated 
at the front of the propagating NZ is also higher. As a 
result, other things being equal, over the same period of 
time in adiabatic conditions the peak temperature rises 
higher in the stack of conductors with high contact resis- 
tance. 

However, it is important to emphasize that neither the 
speed of NZP alone, nor the rate of increase of the peak 
temperature alone can be considered as the good mea- 
sures of the quench protection quality of the conductors. 
If the detection of the quench is based on the voltage that 
builds up across the coil, the higher propagation speed 
of NZ leads to earlier detection, so that the peak tem- 
perature by the time of detection may in fact be lower 
in the stack made out of coated conductors with high 
contact resistance. The rate at which the peak tempera- 
ture increases with increasing voltage drop dT pea k/dV is 
the relevant measure of the quench protection quality of 
the conductors. It is desirable to make this rate as small 
as possible, because it leaves more time for the quench 
protection system to be engaged without the danger of 
overheating the conductor. 

Figure 3 presents the quantitative results obtained by 
solving Eqs. (20) and (21). In Fig. 3(a) the curves 
enumerated (1), (2), and (3) show how the voltage drop 
across the coil, Eqs. (24) and (25), changes with time. 
The curve (1) corresponds to the lowest contact resis- 
tance (see Eq. (23)) A 2 /? 2 . = °- 1 - The value of the 
contact resistance R = 5.1 x 10~ 6 f2 cm 2 corresponds to 
Rq = 5.1 x 10~ 5 fi cm 2 (Eq. (35)) and it is given here for 



illustration purpose only. To the best of our knowledge, 
the value of the contact resistance in currently manu- 
factured coated conductors at T = 4.2 K has not been 
measured. At T = 77 K the contact resistance is much 
smaller than Rq [12j. The curves (2) and (3) correspond 
to progressively greater contact resistance. The scale on 
the right and the upper scale correspond to Vq = 2 mV 
and 7 w 33 a -1 , Eqs. (36) and (31). 

An approximately quadratic time dependence of V(t) 
indicates that over the time interval shown in the Fig. 3 
the NZ expands simultaneously in the lateral and trans- 
verse directions. It is evident that the initiation of the 
NZ growth in the outermost conductors of the stack takes 
place due to the transverse heat flux. The same conclu- 
sion was drawn in Ref. [5|. The characteristic steps in 
dV(i)/dt, most evident in curve (3), correspond to the 
time instances when the two new outermost conductors 
(symmetrically located relative to the center of NZ) reach 
temperature 2\ and start generating voltage. 

Figure 3(b) shows the dimensionless peak temperature 
inside the stack (at the point of NZ origination, Fig. 2(a)) 
as a function of time. For illustration purposes the scale 
on the right shows the peak temperature in absolute units 
for the set of material and operational parameters chosen 
in the text as an example 



T, 



peak 



= Ti 



'peak 



(T c - 7\); T 1 = T C -T 1 = 47 K. (40) 



It is important to note that the temperature difference 
between the curves (2) and (3) is relatively small even 
though curve (3) corresponds to four times greater con- 
tact resistance. In the conductor like the one shown in 
Fig. 1 the excess of heat is generated only at the front of 
the propagating NZ where the current exchange between 
the superconductor and stabilizer takes place across the 
interface between them. Behind the front, the current 
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FIG. 3: The results of the simulation of NZP for three differ- 
ent values of the contact resistance. Curve (1) corresponds to 
the lowest contact resistance, curve (3) to the highest, (a) The 
voltage drop (Eq. (24)) vs time. The scale on the right and 
the upper scale correspond to Vo = 2 mV and 7 = 33 cm~ 
respectively, (b) The peak temperature in the center of NZ 
vs time. The scale on the right corresponds to the set of op- 
erating parameters shown in Eq. (40). (c) A combination of 
the data shown previously - the peak temperature vs voltage 
drop. The dashed line illustrates that when the voltage drop 
reaches a certain threshold, at which point the quench pro- 
tection system may be engaged, the peak temperature is the 
lowest in the stack made out of conductors with the highest 
contact resistance. 



flows only through the copper stabilizer generating the 
minimum amount of heat determined by the value of cur- 
rent and the stabilizer resistance. The peak temperature 
of the initial condition has to be high enough to trigger 
NZ propagation. As the result, the initial phase of the 
time evolution is characterized by the falling temperature 
in the center of the stack due to the energy outflow from 
the point of NZ initiation. 

Combining the results shown in the Figs. 3(a,b) we 
obtained the dependence of the peak temperature in the 
stack on the potential drop across the coil. Here the 
benefits of the increased contact resistance become obvi- 
ous. The slowest temperature rise with increasing voltage 
drop is attained in the stack composed of the conduc- 
tors with the highest contact resistance. As an illustra- 
tion, the vertical dashed line corresponds to a presumed 
threshold which triggers the quench protection system 
(V/V = 100; V = 200 mV). At this level of potential 
drop the peak temperature inside the stack is substan- 
tially lower for the conductors with the highest contact 
resistance because the threshold value is reached much 
sooner. 



B. Insulation with high heat conductance 

(k± = 12). 

The second set of data shown in Fig. 4 corresponds 
to a substantially higher coupling constant k± = 12 (Eq. 
(34)). The large coupling does not necessarily mean some 
anomalously high thermal conductivity of the insulation. 
It can be a result of low resistance of the stabilizer. Since 
the resistivity of copper strongly changes with tempera- 
ture, the real thermal coupling between the turns of a coil 
changes somewhere within the range 0.2 < k± < 12. In 
Fig. 4 all notations and the values of contact resistance 
are the same as in Fig. 3. The only difference is the value 
of the thermal coupling between the neighboring conduc- 
tors in the stack. The higher thermal coupling leads to a 
higher speed of NZ propagation in the transverse direc- 
tion. The normal zone cuts very rapidly through all 21 
conductors in the stack that were considered in the sim- 
ulation. After that, the growth of NZ continues only in 
the lateral direction with a constant speed, which is why 
the time dependence V(t) in Fig. 4(a)is mostly linear, as 
opposed to quadratic in Fig. 3(a). 

The temperature of the initial perturbation has to be 
much higher than that in the case of kj_ — 0.2 in order 
to trigger the NZ propagation. This is why in Fig. 4(b) 
the initial temperature drop is more pronounced than 
that in Fig. 3(b). Figure 4(c) shows the same overall 
result as Fig. 3(c) . The peak temperature in the stack 
made out of conductors with the highest contact resis- 
tance is the lowest for a given value of the voltage drop. 
The dashed line illustrates the peak temperature for a 
presumed threshold V/V = 300 {V = 600 mV). 
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FIG. 4: The results of the simulation of NZP for three differ- 
ent values of the contact resistance and the thermal coupling 
ki = 12. The curve (1) corresponds to the lowest contact 
resistance, the curve (3) to the highest, (a) The voltage drop 
(Eq. (24)) vs time. The scale on the right and the upper scale 
correspond to Vb = 2 mV and 7 = 33 cm -1 respectively, (b) 
The peak temperature in the center of NZ vs time. The scale 
on the right corresponds to the set of operating parameters 
shown in Eq. (40). (c) A combination of the data shown pre- 
viously - the peak temperature vs voltage drop. The dashed 
line illustrates that when the voltage drop reaches a certain 
threshold, at which point the quench protection system may 
be engaged, the peak temperature is the lowest in the stack 
made out of conductors with the highest contact resistance 



V. SUMMARY AND RUMINATION OVER 
QUENCH PROTECTION 

j previous 
quench in coated conductors [0, 0| to a case of a stack 
of conductors imitating a part of a pancake coil. Our 
criterion of the quench protection quality of a conductor 
is the rate at which the peak temperature inside the nor- 
mal zone increases with the voltage drop across the coil. 
The smaller this rate is, the more reliable the quench 
protection system can be made. The detection of the 
quench is based on the voltage monitoring. To avoid 
the false alarms the threshold of triggering the quench 
protection system has to be made high enough. There- 
fore, it is highly desirable that at the moment when the 
threshold potential is reached the peak temperature in- 
side NZ was as low as possible, certainly well below the 
temperature above which the conductor may suffer an 
irreversible damage. The results shown in Figs. 3 and 
4 clearly indicate that the increased contact resistance 
leads to a faster build up of the voltage drop which more 
than compensates for the increased heat dissipation at 
the front of the propagating NZ. 

The origin of the increased longitudinal speed of NZP 
(along the conductor) was explained in 7]. The cur- 
rent exchange length A, Eq. (17), replaces the thermal 
diffusion length as the dominant length scale that deter- 
mines the NZP speed. In a coil or a stack like the one 
shown in Fig. 1(b) the transverse rate of NZP also in- 
creases due to increased contact resistance. The results 
of our simulation, as well as the experiments and simula- 
tions presented in Ref. [5| show that the propagation of 
the NZ in the transverse direction is caused by the heat 
transfer across the insulation, rather than by the NZP 
along the conductor. The coated conductors with high 
contact resistance are less stable [7j, so that it takes less 
energy (and shorter time) to trigger NZ in the outlying 
sections of the coil. Therefore, in a coil the advantages 
of the increased contact resistance in increasing the rate 
of potential build up are even more pronounced than in 
a linear conductor cooled from the surface. This is due 
to a symbiotic effect of increased longitudinal NZP speed 
and reduced stability. 

The main apparent shortcoming of the suggested ap- 
proach to improving the quench protection quality of 
coated conductors is that by increasing the contact resis- 
tance we artificially reduce the stability margins of the 
device with respect to quenchQ- On the surface, the 
high stability of coated conductors is a plus. It makes 
quenches less frequent. However, no magnet designer will 
be so bold as to suggest that a coil made out of coated 
conductors does not need a quench protection system. 
The prospect of a Black Swan quench like a recent one 
at the Large Hadron Collider [13| (unrelated to coated 
conductors, but noteworthy in the context of this paper) 
makes it inevitable that no matter how stable the current 
generation of coated conductors appears, a full-fledged 
protection system has to be designed and built. Then, 



what is the advantage of the high stability conductor? 
None - in terms of the cost of the design and construc- 
tion of the quench protection system. Moreover, because 
of the low NZP speed in the current generation of coated 
conductors such a system will require a greater effort to 
design and validate than the state-of-the-art quench pro- 
tection systems in existence today. The low-T c supercon- 
ducting magnets frequently quench when the current is 
ramped up for the first time. During these training ses- 
sions the weak links in the coil are revealed and healed. 
A very stable magnet might not quench at all for the first 
time it is energized with the result that a weak spot will 
remain undetected until it reveals itself catastrophically 
at some point during the lifetime of the device. 

It is obvious then, that a way to overcome these prob- 



lems is to reduce the stability margins of the conductor by 
inserting a very thin (thermally transparent), but high re- 
sistance interface between the superconducting film and 
the stabilizer. Reduced stability and increased speed of 
NZP in such a conductor will revert the potential Black 
Swan quenches into already familiar, more frequent, but 
tamed variety. It is possible, of course, that the exist- 
ing systems of quench detection may be upgraded and 
refined to a level that will be sufficient for use with the 
current generation of coated conductors. However, if it 
turns out not to be the case, then this and the previous 
publications jy, 0] may be considered as a rough draft of 
the proverbial Plan B. 
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